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ABSTRACT 


An  equilibrium  model  of  a manpower  system  is 
developed  based  on  the  notion  of  a career  flow. 
Institutional  constraints  and  measures  of  system 
performance  are  linear  functions  of  the  career 
flow.  A typical  optimal  design  problem  is  formu- 
lated and  a solution  procedure  is  developed.  The 
optimization  problem  is  a generalized  linear  pro- 
gram in  which  columns  are  generated  by  solving  a 
shortest  path  problem.  Upper  and  lower  bounds  on 
the  optimal  value  function  can  be  developed  at 


OPTIMAL  DESIGN  OF  A MANPOWER  SYSTEM 


0.  INTRODUCTION 

This  paper  presents  a model  for  the  design  of  a manpower  system  and 
develops  a procedure  for  Its  solution.  The  manpower  system  Is  assumed  to 
be  In  temporal  equilibrium  and  to  be  a strict  hierarchy;  l.e.,  a person 
In  rank  1 either  stays  In  rank  1 , moves  to  rank  1 + 1 , or  leaves  the 
system.  The  equilibrium  assumption  Is  critical  to  the  model.  The  assump- 
tion of  strict  hierarchy  Is  made  for  simplicity;  a more  complicated  system 
could  be  considered  In  the  same  way  with  more  complicated  results. 

The  approach  to  modeling  the  manpower  system  was  suggested  by  the  equil- 
ibrium chain-flow  model  of  Oliver  and  Hopkins  [7].  Their  Idea  of  chain  flow 
Is  exactly  the  notion  of  career  flow  used  In  this  paper.  Oliver  and  Hopkins 
stipulate  a set  of  possible  student  and  faculty  flow  paths  through  a univer- 
sity and  then  calculate  the  path  flows  that  are  consistent  with  meeting  in- 
stitutional constraints  at  minimum  cost.  This  paper  presents  these  Ideas  In 
a more  general  setting  and  shows  how  the  problem  of  having  a very  large  num- 
ber of  possible  career  paths  can  be  handled. 

The  motivation  for  the  procedure  developed  in  this  paper  can  be  found 
in  Schmidt,  et  al.  [8],  [9].  The  problem  posed  In  [8]  and  [9]  is  to  design 
a manpower  system  In  order  to  maximize  the  effectiveness  to  cost  ratio  sub- 
ject to  several  institutional  constraints.  A rather  complicated  nonlinear 
optimization  approach  that  deals  with  an  aggregate  version  of  the  problem 
is  presented  In  [8].  In  contrast  the  approach  described  here  Is  flexible, 
general,  deals  with  the  entire  problem,  and  Is  easily  solved. 

The  problem  of  system  design  should  be  contrasted  with  the  problem  of 


system  control.  In  [4],  [5],  procedures  for  controlling  graded  manpower 
systems  were  developed.  These  procedures  assumed  that  certain  ratios  among 
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manpower  stocks  and  flows  were  known,  stationary,  and  governed  the  manpower 
flow  process.  In  contrast,  the  approach  developed  in  this  paper  actually 
designs  the  flow  process  and  yields  the  ratios  as  policy  variables. 

Section  I describes  the  manpower  system  and  the  equilibrium  flow  equa- 
tions. The  notion  of  a career  is  Introduced  in  Section  II  and  several  im- 
portant characteristics  of  the  manpower  system  are  expressed  in  terms  of 
career  flows.  In  Section  III  linear  programming  optimization  model  is  pre- 
sented and  Section  IV  describes  a column  generation  technique  for  its  solu- 
tion; the  columns  are  generated  by  solving  a shortest  path  problem.  In 
Section  V a termination  criterion  is  developed.  The  special  structure  of 
the  problem  allows  one  to  determine  if  the  current  solution  is  within  a 
specified  percentage  error  of  the  optimal  solution.  A subsequent  paper  will 
deal  with  applications  of  the  techniques  developed  here  to  the  enlisted  force 
of  the  U.S.  Navy. 
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I.  THE  SYSTEM 

The  manpower  system  is  an  N-rank  hierarchy  that  is  observed  at  equally 
spaced  points  in  time  t « 0,1,2,  ...  . The  time  Interval  (t  - l,t]  , after 
time  t - 1 and  up  to  and  including  time  t , is  called  period  t . At  ob- 
servation time  t each  individual  is  given  a classification  (i,u)  , accord- 
ing to  their  rank,  1 , and  period  of  service,  u . A period  of  service 
equal  to  u at  time  t indicates  the  individual  entered  the  system  in  period 
t + 1 - u : l.e.  in  the  time  Interval  (t  - u,t  + 1 - ul  . An  individual's 
length  of  service  is  the  number  of  aorrrpleted  periods  of  service;  thus  an  in- 
dividual in  class  (l,u)  has  length  of  service  u - 1 . The  maximum  length 
of  service  is  M - 1 ; thus  an  individual  in  class  (i,M)  must  leave  the 
system. 

Figure  1 shows  the  possible  movements  from  class  (l,u)  . The  individual 
in  (l,u)  moves  either  out  of  the  system  (to  class  0),  stays  in  rank  1 , or, 
if  1 < N , is  promoted.  There  are  no  demotions  or  double  promotions,  and 
each  person  enters  in  class  (1,1)  . 


FIGURE  1:  Possible  Transitions  from  Class  (l,u). 
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We  assume  the  system  Is  In  equilibrium.  Let  s(l,u)  be  the  number 
of  people  In  class  (l,u)  , and  let  f(l,u,j)  be  the  number  of  people  per 
period  who  move  from  class  (l,u)  to  class  (j,u  + 1)  . Notice  the  assump- 
tion on  possible  movements  Implies  f(l,u,j)  = 0 If  j Is  not  equal  to 
0 , 1 , or  1 + 1 . Of  course,  f(N,u,N  + 1)  ■ 0 , since  there  are  only  N 
ranks,  and  f(l,M,l)  ■ f(l,M,l  + 1)  = 0 since  one  must  leave  the  system 
during  the  period  of  service. 

The  equilibrium  flow  equations  are 


(2)  (1)  for  u = 1 , 1 = 1 

s(l,l)  = f(l.l,0)  + f(l,l,l)  + f(l,l,2)  . 

(11)  for  l<u<M,  1-1 

f(l,u  - 1,1)  - s(l,u)  - f(l,u,0)  + f(l,u,l)  + f(l,u,2)  . 

(Ill)  for  l<u<M,l<l<N 

f(l  - l,u  - 1,1)  + f(l,u  - 1,1)  =•  s(l,u)  » f(l,u,0)  + f(l,u,l) 

+ f(l,u,u  + 1)  . 

(Iv)  for  l<u<M,  1 = N 

f(N  - l,u  - 1,N)  + f(N,u  - 1,N)  - s(N,u)  “ f(N,u,0)  + f(N,u,N)  . 
(v)  for  u - M , 1-1 

f(l,M  - 1,1)  - s(l,M)  - f(l,M,0)  . 

(vl)  for  u ■ M , 1 < 1 

f(l  - 1,M  - 1,1)  + f(l,M  - 1,1)  - s(l,M)  - f(l,M,0)  . 


The  equilibrium  flow  for  N - 3 , M - 5 Is  shown  In  Figure  2,  retirement 
flows  (l,u,0)  are  not  shown. 

The  variable  s(l,u)  and  f(l,u,J)  allow  for  a measure  of  the  effec- 
tiveness and  cost  of  different  equilibrium  flow  patterns.  Let  e(l>u)  be 
the  effectiveness  of  an  Individual  In  class  (l,u)  , c(l,u)  the  cost, 
p(l,u)  the  cost  associated  with  a promotion  from  (l,u)  to  (1  + l,u  + 1) 
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and  r(i,u)  the  cost  of  a retirement  from  (l,u)  . The  cost  and  effective- 
ness are  respectively: 

(3)  (i)  I } c(i,u)s(i,u)  + p(i,u)f (i,u,i  + 1)  + r(i,u)f (i,u,0)  . 

i u 


I 1 e(i,u)s(i,u) 
1 u 


Although  Equations  (2)  and  (3)  describe  the  equilibrium  conditions  and 
possible  criteria  they  are  not  in  a suitable  form  to  describe  the  structure 
of  allowable  flow  patterns.  That  structure  will  be  discussed  in  the  next 
section. 
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II.  CAREERS 

A career  is  a path  through  the  manpower  system  that  commences  at  rank 
1 with  u = 1 . Let  k index  the  career.  Each  career  has  several  im- 
portant characteristics.  Some  of  these  are: 

(4)  (i)  n(i,k)  ~ The  number  of  times  counted  in  class  i . 

(ii)  m(i,k)  ~ Attainment  of  class  i : this  variable  is  one  if  career 

k includes  class  i (i.e.  if  n(i,k)  > 0),  zero  if 

career  k does  not  include  class  i : i.e,  n(i,k)  * 0 

(iii)  a(i,k)  ~ Length  of  service  upon  entry  to  class  i . 

(iv)  d(v,k)~  The  number  of  periods  with  length  of  service  greater  than 

V ; i.e.  time  spent  in  the  professional  force. 

(v)  e(k)  ~ The  total  effectiveness  of  career  k . 

(vl)  c(k)  ~ The  total  cost  of  career  k . 

(vii)  n(k)  ~ Total  length  of  career  k , n(k)  = (n(i,k)  . 

i 

A career  may  be  characterized  by  transitions  between  nodes.  Let 

. ( 1 if  career  k moves  from  (l,u)  to  (j,u  + 1) 

rU.u.j)  - 

' 0 otherwise . 

The  characteristics  described  in  (4)  can  be  expressed  in  terms  of  the 
variables  f^  . 

(5)  (i)  n(i,k)  - I I f'^d.u.j)  . 

u j 

(11)  m(l,k)  - I f*^(l  - l,u,l)  for  1 > 2 . 
u 

(111)  a(l,k)  - I uf'^d  - l,u,l)  for  1 > 2 . 
u 

(lv)d(v,k)  - I I I f^(l,u,j)  . 
u>v  1 j 
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(v)  c(k)  = I I I c(i,u, j)  . 

u i J 

(vl)  e(k)  = I I e(i,u)  J . 

u i j 

(vii)  n(k)  “III  . 

i u j 


The  costs  c(i,j,j)  are  calculated  from  basic  costs  described  in  Section  I. 


(6)  (i)  c(i,u,0)  = c(i,u)  + r(i,u)  . 

(ii)  c(i,u,i)  = c(i,u)  . 

(iii)  c(i,u,i  + 1)  = c(i,u)  + p(i,u)  . 


The  special  structure  of  the  organization  implies  that  f (i,u,i  + 1)  = 1 
for  at  most  one  value  of  u . Thus  at  most  one  term  in  the  summations  (5:ii) 
and  (5: iii)  is  equal  to  one. 

Let  g(k)  be  the  number  of  people  starting  career  k each  period. 

Under  the  steady-state  assumptions  there  will  be  g(k)  individuals  in  each 
phase  of  the  career.  It  follows  that; 


(7)  (i)  I n(i,k)g(k)  is  the  total  number  of  people  in  class  i . 

k 

(ii)  I m(i,k)g(k)  is  the  number  of  entrants  that  ever  reach  class  i . 
k 

(iii)  I a(i,k)g(k)  is  the  total  length  of  service  of  individuals 
k 


entering  class  i . 

(iv)  J]d(v,k)g(k)  is  the  number  of  people  with  length  of  service 
k 

greater  than  v . 

(v)  I c(k)g(k)  is  the  total  cost  incurred, 
k 

(vi)  [ e(k)g(k)  is  the  total  effectiveness, 
k 

(vll)  I n(e)g(k)  is  the  total  number  of  people  in  the  system, 
k 


This  description  of  the  system  is  based  on  the  notion  of  a career.  The  next 
section  demonstrates  how  the  system  can  be  designed  in  an  optimal  manner. 
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III.  THE  OPTIMIZATION  PROBLEM 

The  model  can  be  used  with  several  objectives.  It  is  possible  to  maxi- 
mize effectiveness  subject  to  a cost  constraint.  An  alternative  is  to  mini- 
mize cost  subject  to  an  effectiveness  constraint.  A third  objective  is  to 
maximize  an  effectiveness/cost  ratio.  These  objectives  are  discussed  at 
length  in  an  appendix.  The  choice  of  an  objective  does  not  change  the  op- 
timization procedure.  The  example  presented  below  is  to  minimize  cost  sub- 
ject to  effectiveness  and  system  constraints. 

Let  q(i,k)  = a(i,k)  - a(l)m(l,k)  be  excess  experience  upon  entry  to 
rank  i . The  system  constraints  are 


(8)  (i) 

(li) 

(ili) 

(iv) 

(v) 

(vl) 


J n(i,k)g(k)  > r(i)  , at  least  r(i)  people  in  rank  1 ; 
k 

i = 1,2 N . 

^ q(i,k)g(k)  > 0 , the  average  length  of  service  upon  entry  to 
k 

class  i at  least  a(i)  for  i = 2,  ...»  N . 

^ d(v,k)g(k)  > d , a minimum  size  of  the  professional  staff, 
k 

'l  e(k)g(k)  > e a minimum  level  of  effectiveness, 
k 

^ n(k)g(k)  = r ; fixed  total  strength, 
k 

^ c(k)g(k)  the  cost, 
k 


The  linear  program  is  to  minimize  (8:vl)  subject  to  the  constraints  (8;l-v). 
This  linear  program  has  2(1  + N)  constraints  and  a nonnegative  variable, 
g(k)  for  each  career.  The  next  section  describes  a method  for  solving  this 
problem. 

Equation  (8)  describes  several  examples  of  constraints.  There  are  other 
possibilities:  three  additional  examples  are  presented  below.  The  first  is 
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r 


5i{m(i  + l,k)  - a(l  + l,i)m(i,k)}g(k)  > 0 . 
k 


The  fraction  of  those  who  reach  rank  i + 1 given  they  have  reached  rank  i 
exceeds  a(l  + 1,1)  . A second  example  is 


j;{m(i  + l,k)  - 6(1  + l,i)n(i,k)}g(k)  > 0 . 
k 


The  fraction  of  those  in  rank  i that  are  promoted  to  rank  1+1  per 
year  is  at  least  6(1  + l,i)  . Finally, 


)]{d(h,k)  > ?d(Jl,k)}  > 0 . 
k 


The  fraction  of  those  with  length  of  service  greater  than  k who  remain  to 
have  length  of  service  greater  than  h > is  at  least  ^ . 


>4 
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IV.  SOLUTION  PROCEDURE 

The  solution  procedure  is  outlined  in  this  section.  It  is  a column 
generation  algorithm:  see  [1],  [2],  and  [3]  for  the  origins  and  applica- 

tions of  this  idea.  Column  generation  is  a variant  of  the  revised  simplex 
method.  The  columns  of  the  activity  matrix  are  not  stored  explicitly. 

They  are  generated  by  solving  a sub-problem. 

The  section  describes  one  iteration  of  the  procedure.  Assume  a feasible 
basis  K is  available.  Let  k index  the  career  in  the  bases.  There  are 
simplex  multipliers  associated  with  this  basis;  each  multiplier  is  associated 
with  a constraint. 

(9)  (i)  ’l'(i)  ~ stock  constraints  i = 1.  ...,  N . 

(il)  A(i)  ~ average  length  of  service  constraints  i = 2,  ...»  n . 

(iii)  p ~ professional  staff  constraint. 

(iv)  Y ~ effectiveness  constraint. 

(v)  a ~ total  size  constraint. 


These  simplex  multipliers  satisfy  the  following  linear  equality  for  careers 
k in  the  current  bases. 

(10)  ^ 4'(l)n(i,k)  + ^ A(l){q(i,k)>  + d(v,k)  + Ye(k)  + an(k)  + c(k)  = 0 . 

i i>2 

The  current  basis  is  optimal  if  (i|^fX,y,Y)  are  nonpositive  and  for  all  careers  k , 

(11)  ^ 4'(i)n(i,k)  + ^ X(i)q(i,k)  + d(v,k)  + Yc(k)  + an(k)  + c(k)  > 0 . 

i i>2 

a 
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From  (10)  note  that  the  optimal  value  of  the  subproblem  (12)  is  less 
than  or  equal  to  zero.  If  the  optimal  value  of  the  subproblem  is  zero,  then 
(11)  is  satisfied  and  the  current  basis  is  optimal.  If  the  optimal  value  of 
the  subproblem  is  negative,  then  a new  career  can  be  found  that  is  eligible 
to  enter  the  current  basis.  The  procedure  for  solving  (12)  and  generating 
the  new  career  is  presented  below. 

Substitute  the  Equations  (5),  from  Section  2,  into  (12)  and  combine 
terms.  This  yields 

(13)  Minfn  I l(i,u,J)f^(l,u,j)l  . 

k [u  j 1 J 

The  correct  values  of  1(1, u,j)  are  given  below  where  d(u)  ■ 1 if  u ^ v , 
otherwise. 

(lA)  (i)  f(l,u,l)  =>  c(l,u)  + I'd)  + iJd(u)  + ye(i,u)  + a . 

(ii)  4(i,u,0)  ■ i(i,u,i)  + r(i,u)  . 

(ill)  a(i,u,l  + 1)  - i(i,u,i)  + p(l,u)  + X(i  + 1) (u  - a,i  + 1) 
for  i < N - 1 . 

a 

Problem  (13)  is  a shortest  path  problem.  The  network  was  described 
in  Figure  1,  Section  1.  The  l(i,u,j)  are  the  lengths.  Let  V(l,u)  ■ 
length  of  the  shortest  path  starting  from  node  (l,u)  out  of  the  system. 

v(i,M)  - a(i,M,o)  . 

For  u < M 

!Ui,u,0) 

Jl(i,u,l)  + v(i,u  +1)  If  i < N 

i(i,u,l  + 1)  + v(i  + l,u  + 1)  . 

( il(N,u,0) 

(ii)  v(N,u)  - Min  { 

( )l(N,u,N)  + v(N,u  + 1)  . 
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Finally  v = Min  [v(i,l)]  is  the  length  of  the  shortest  path.  If  v*  < 0 , 
the  shortest  path  defines  a career  that  is  eligible  to  enter  the  bases.  If 

A 

V *=  0 , the  current  solution  is  optimal. 
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V.  BOUNDS  ON  THE  OBJECTIVE 

The  column  generation  scheme  described  In  Section  IV  can  be  modified 
slightly  to  give  a series  of  upper  and  lower  bounds  on  the  optimal  value 
of  the  objective.  Once  a feasible  solution  has  been  obtained  the  cost  of 
the  current  feasible  solution  Is  an  upper  bound  on  the  optimal  value.  Let 
U be  this  upper  bound. 

N 

(16)  U - - \ i/)(l)r(i)  + pp  + ye  + or  . 

1=1 

Recall  that  the  shortest  path  calculation  solves  the  problem 

(17)  Min  Fc(k)  + \ (^(i)n(i,k)  + \ X(l)q(l,k)  + pp(k)  + ye(k)  + on(k)'l  . 

k L 1 1>2  J 

Without  loss  of  generality,  'll  , X , v , y are  nonpositive.  If  they  were 
positive  a slack  variable  would  have  been  admitted  to  the  basis.  Let 
<t>  “ ('('»X,p,y,a)  , and  call  the  minlmand  in  (17),  x(k,i|i)  . The  lower 
bound  is  developed  by  studying  a fractional  program  (18).  Notice  that 


(18) 


the  fractional  optimization  problem  (18),  will  have  an  optimal  value  of  the 
same  sign  as  (17).  Let  g (k)  be  the  number  of  people  following  career  k 
in  an  optimal  solution  of  the  overall  problem  and  let  V be  the  optimal 
value  of  the  overall  problem. 

V = I c(k)g*(k)  < U . 
k 


(19) 
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Now  define  z(k)  = {n(k)g  (k)}/r  ; z is  a nonnegative  vector  that  sums 
to  one.  Suppose  career  h solves  (18),  with  value  W , then  for  all  k 


(20) 


x(h,(<.)  , x(k,()>) 
n(h)  - n(k) 


Multiply  (20)  by  z(k)  and  sum;  this  yields 


(21) 


W 5 7 I c(k)g  (k)  + I I//.  I n(i,k)g*(k)  + I X(l)  I q(l,k)g*(h)  + 
" 1_  i ^ k i>2 


U I p(k)g  (k)  + Y I e(k)g  (k)  + a ^ n(k)g  (n) | . 
k k k 


) 


From  the  feasibility  of  g and  the  nonpositivity  of  ('I'»^»M»Y)  • this 
becomes 


(22) 


W < 7 [V  - U]  or  U + rW  < V . 


If  e is  a desired  percentage  error  in  the  calculation  of  V , then  the 
stopping  rule  is  to  terminate  if 


(23) 


< c 


Notice  that  calculation  of  this  lower  bound  involves  the  solution  of 
(18)  rather  than  (17).  However,  (17)  is  a rather  easy  to  solve  shortest 
path  problem  and  (18)  is  a more  difficult  to  solve  fractional  program.  It 
is  possible  to  approximate  the  '.falue  of  W using  (17).  Let  Y(k,i(i,u)) 
and  L(u)  be  defined  as 

Y(k,i^,u))  ■ x(k,4>)  - u)n(k) 

L((o)  - Min  [Y(k,4.,a))l  . 
k 


(24) 


16 


F 


Notice  that  evaluating  L((d)  is  equivalent  to  solving  problem  (17)  with 
a replaced  by  (a  - w)  . The  optimal  solution  of  (18)  is  found  by  finding 

* 4t  * 

the  0)  that  satisfies  L(o)  ) = 0 , in  fact  u » W . The  reader  can  show 

that  L((jj)  is  a decreasing,  piecewise  linear,  concave  function  of  u . 

Furthermore,  suppose  h is  the  career  that  solves  (17).  If  there  is  no 

unique  optimal  career,  then  h is  selected  to  minimize  n(h)  among  all 

shortest  paths.  Problem  (17)  evaluates  L(0)  , and  the  left  derivative  of 

L at  w = 0 , see  [3],  is  -n(h)  . Therefore  = L(0)/n(h)  is  an  upper 
* 

bound  for  u . Define  u>2  by  = -cU/r  ; notice  from  (23)  and  the  mono- 

* 

tonicity  of  L that  cutoff  occurs  if  and  only  if  ooj  < m 

The  three  possible  cases  that  can  occur  are  shown  in  Figure  3.  First, 
in  3. a,  0)  cutoff  is  possible.  If  0)2  ^ , then  L(u2) 

is  evaluated  by  solving  (17)  with  a replaced  by  a - . In  Figure  3b., 

L(u2)  < 0 , therefore  m and  cutoff  does  not  occur.  In  Figure  3c. 

it 

L(to_)  > 0 , so  u-  < u)  and  calculations  are  terminated  with  a percentage 

2.2=  ^ ^ 

error  less  than  e . ' 


i 


i 
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